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1. Motivation 

Even though the Green's (correlation) functions of pure Yang-Mills theories are not per se 
physical objects (for they show an explicit dependence on the gauge-fixing and renormalization 
scheme used), the understanding (or lack thereof) of their infra-red (IR) properties has become an 
increasingly interesting topic in the last five years. The interest in these studies lies in the fact that 
reliable information on their non-perturbative structure is essential for unraveling the IR dynamics 
of QCD in general, and the mechanism behind confinement in particular. 

The exploration of the IR sector of such theories is currently pursued through mainly two non 
perturbative tools, namely the lattice - where space-time is discretized and the quantities of inter- 
est are evaluated numerically - and the Schwinger-Dyson equations (SDEs) - the infinite set of 
integral equation governing the dynamics of the Green's functions. Each of the aforementioned 
methods has its own weaknesses, some of these being common to both (e.g., the need to fix a 
gauge, which in the lattice calculations generally coincides with the Landau gauge), some other 
being instead specific to the method under scrutiny. In the lattice case, for instance, one has to keep 
under control the various sources of systematic errors (e.g., discretization effects, finite volume ef- 
fects, Gribov copies effects), while at the same time providing enough computational power (that is 
large volume lattices) to study the deep IR region; in the SDEs case, the problem lies in devising a 
self-consistent truncation that would reduce the infinite tower of equations to a manageable subset, 
without introducing artifacts distorting the properties one has endeavored to investigate. Recently, 
a lot of progress has been made in addressing several of the aforementioned problems in both meth- 
ods [0» ^ e||> and we definitely have reached a point where SDE predictions can be systematically 
compared against lattice results. 

Over the last few years, in fact, high quality ab-initio lattice gauge theory computations have 
established beyond any reasonable doubt that the gluon propagator and the ghost dressing function 
of pure Yang-Mills in the Landau gauge saturates in the deep IR at a finite, non-vanishing value. 
The emergence of these massive solutions from the lattice together with their confirmation at the 
level of SDEs [Q, ^j, has caused a paradigmatic shift in our current understanding of QCD; in 
particular, the gluon mass generation scenario, put forth by Cornwall and others in the eighties [§], 
appears to be the preferred interpretation of what the underlying IR dynamics might be 1 . In this 
latter picture the gluon acquires dynamically an effective (momentum-dependent) mass, which 
accounts for the IR finiteness of the gluon propagator. The local gauge invariance and (after gauge 
fixing) the BRST symmetry remain intact, because the original QCD Lagrangian is not altered 
at any point; indeed, the non-perturbative mechanism responsible for the mass generation is the 
four-dimensional generalization of the celebrated Schwinger mechanism p0|]. 

In Fig. [l| we show a compilation of lattice data for the Landau gauge SU(2) [[I], [IT]] and 
SU(3) Jl2| ] gluon propagator (top panels) and ghost dressing function (bottom panels) from two 
independent groups. Although for each group the various simulation parameters vary (including 
lattice spacing, gauge group employed, and gauge fixing algorithm), all results shares the common 

'The SDEs admit a second class of solutions that go under the name of scaling solutions. They are characterized 
by an IR vanishing gluon propagator and, correspondingly, an IR diverging ghost dressing function |Q, and have been 
predicted in confinement scenarios such as the one of Kugo-Ojima ^ and Gribov-Zwanziger |^]. So far, however, these 
solutions have not been found in (Landau gauge) lattice simulations. 
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feature of the appearance of a plateau for the gluon propagator in the deep IR region, which is one 
of the most salient and distinctive predictions of the the gluon mass generation mechanism to be 
discussed below. In addition, the ghost dressing function show no enhancement in the deep IR; 
instead, it again saturates to a constant. 

Now, define the full gluon and ghost propagators in the gauges as 



A M v(?) 



^ v {q)=P^(q)Kq 2 ), D{q 2 )=i 



F(q 2 



(1.1) 



where P^ v (q) = gfiv(q) — q^qv/q 2 is the transverse projector, and the scalar function A(q 2 ) is re- 
lated to the gluon self-energy H^ v (q) = P flv (q)Tl(q 2 ) by A~ l (q 2 ) = q 2 + iH(q 2 ); F(q 2 ) above is 
the ghost dressing function, which is related to the ghost self-energy £(q 2 ) through q 2 F~ 1 (q 2 ) = 



q" — \t{q r). Given the above, one observes that the (Landau gauge) lattice data sets can be accu- 
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where Ca is the Casimir eigenvalue of the gauge group's adjoint representation [Ca = N for SU(A0], 
while p, pi 2 and g\ t 2 are treated as free fitting parameters 2 . From the expressions above we clearly 
see the anticipated dependence on the momentum transfer q 2 displayed by the dynamically gener- 
ated mass; in particular, the latter assumes a non-zero value mo in the IR while dropping rapidly in 



the ultra-violet (UV) in a way consistent with the operator product expansion results Q14Q (see next 
section). Clearly mo acts as a physical mass scale: it regulates the perturbative renormalization 
group logarithm, so that, instead of diverging at the Landau pole, it saturates at a finite value. For 
the SU(3) case, the best fit is obtained for a gluon mass mo ~ 500 MeV (with the data renormalized 
at jU = 4.3 GeV). 

Lattice studies in the Landau gauge, therefore, rule out the possibility of scaling solutions 
with nontrivial infrared exponents (and consequently the Kugo-Ojima scenario). In addition, the 
original formulation of the Gribov-Zwanziger scenario is clearly disfavored, and must be drastically 
modified through the inclusion of dimension two condensates, in order to be reconciled with the 



lattice data [15] 



In the rest of this talk we will concentrate on the Pinch Technique - Background Field Method 
(PT-BFM) framework [jS |16[ |l7|], where the aforementioned lattice findings may be naturally ac- 
commodated. Indeed, it should be noticed that the discovery of the key underlying ingredient, 
namely the dynamical generation of a gluon mass, coincided historically with the invention of the 
PT long before lattice simulations of QCD Green's functions were even contemplated. 

2. Dynamical gluon mass generation, center vortices and confinement 



As mentioned earlier, already in the fifties Schwinger [jlOj pointed out that the gauge invariance 
of a vector field does not necessarily imply zero mass for the associated particle. Schwinger's idea 
is very simple. Suppose that the vacuum polarization Yl(q 2 ) [where H(q) = q 2 U(q 2 )] acquires, 
for some yet to be specified dynamical reason, a pole in the zero momentum transfer limit q 2 = 0, 
with positive residue \i 2 , i.e., H(q 2 ) = \i 2 jq 2 ■ Then (in Euclidean space) A~ l (q 2 ) = q 2 + jJ. 2 and 
A -1 (0) = jU 2 : that is, the vector meson has become massive, even if the gauge symmetry forbids a 
mass at the level of the fundamental Lagrangian. 

There is no physical principle which would preclude tl(q 2 ) from acquiring such a pole, even 
in the absence of elementary scalar fields; in a strongly-coupled theory like non-perturbative Yang- 
Mills, for example, this may happen for purely dynamical reasons, since strong binding may gener- 



ate zero-mass bound-state excitations Q18Q. The latter act like dynamical Nambu-Goldstone bosons, 
in the sense that they are massless, composite, and longitudinally coupled; but, at the same time, 
they differ from Nambu-Goldstone bosons as far as their origin is concerned: they do not originate 
from the spontaneous breaking of any global symmetry 3 

In the next section we shall see that the way the Schwinger mechanism is incorporated into the 
PT-BFM framework is through a particular form of the vertices appearing in the SDEs, which will 
display dynamical massless poles ~ 1/g 2 - The demonstration of the existence of such a zero-mass 



2 For the SU(3) lattice simulations, fl must be chosen to coincide with the renormalization point, while for the SU(2) 
case it can be treated as an adjustable parameter. 

3 Notice that the usual Higgs mechanism corresponds to a particular case of Schwinger's mechanism, where the 
residue of the pole is saturated by the vacuum expectation value of an elementary scalar field. 
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Figure 2: Self-consistent picture of confinement in the dynamical gluon mass generation scenario. 

bound state, is a difficult dynamical problem, usually studied by means of Bethe-Salpeter equations 
(see, e.g., [|I9|]). In what follows therefore we will rather assume that the theory can indeed generate 
the required poles, and see how far this assumption can lead us. 

First of all, it should be noticed that a (gauge invariant) dynamically generated gluon mass 
cannot be a "hard mass", but will depend non-trivially on the momentum transfer q 2 , and will 
decrease at large momentum (in order to avoid the existence of a corresponding bare mass). Indeed, 



application of the OPE within the PT framework [ |14| ] shows that the gauge invariant contribution 
of the condensate (Tr G^ v ) to the gluon propagator appears in the right way as to be interpreted as 
contributing to the running mass, with the asymptotic ultraviolet behavior m 2 {q 2 ) ~ (Tr G 2 lv )/q 
(actually, powers of logarithms of q can also occur, but we ignore them here). At this point, one 
can also write down an effective low-energy theory describing (gauge invariantly!) the low energy 
dynamics of massive gluons. To be precise, such a theory would be a "quantum" effective theory, 
since it resums all the quantum effects (not present in the classical action) that give rise to the 
dynamical gluon mass. The action we are looking for is the one corresponding to the gauged non- 
linear sigma model and known as "massive gauge-invariant Yang-Mills" p0|], with Lagrangian 
density 

J^mym = \fI v -m 2 Tr [A^-g- l U(9)d fl U-\9)] 2 . (2.1) 

In the formula above = g £ a A a A" the X a are the SU(/V) generators (with TrX a Xb = 2S a b), and 
the. N xN unitary matrix U (6) = exp [iiA fl a l describes the scalar fields 6 a . Note that J2? MYM is 
locally gauge-invariant under the combined gauge transformation 

A'^=VA^V- l -g- l [d^V}V-\ U' = U(e')=VU(9), (2.2) 

for any group matrix V = exp [i^A a ft) a (.x)] , where (O a {x) are the group parameters. 

This action give rise to a plethora of possible "quantum" solitons, i.e., localized finite-energy 
configurations of gauge potentials corresponding to extrema of the effective action: center vortices, 
nexuses and sphalerons. Of these the most important are center vortices, since they show a long- 
range pure gauge term in their potentials, which endows them with a topological quantum number 
corresponding to the center of the gauge group [Z/v for SU(Af)]. This in turn is responsible for 
quark confinement and gluon screening. Specifically, center vortices of thickness ~ , form 
a condensate (Tr G^ v ) ^ because their entropy (per unit size) is larger than their action. This 
condensation furnishes an area law for the Wilson loop in the fundamental representation, thus 



confining quarks Hq, 20]. In addition, the adjoint potential shows a roughly linear regime followed 



by string breaking when the potential energy is about 2mo, corresponding to gluon screening [21]. 
Thus, a dynamically generated gluon mass lead us to the confinement scenario of Fig. || 



5 



Dynamical gluon mass generation and the IR sector of QCD 



Daniele Binosi 




M (05) 



Figure 3: SDE for the conventional gluon self-energy. 

3. Schwinger-Dyson equations in the PT-BFM 

At this point one would like to be more quantitative, and use the SDEs in order to thoroughly 
study the behavior of the gluon propagator and ghost dressing function in the dynamical gluon 
mass generation scenario. The key ingredient enabling such an analysis has been the development 
of a gauge invariant truncation scheme for non- Abelian SDEs [§]. 

To understand where the problem lies, let us consider the conventional SDE for the gluon 
self-energy (Fig. ||). Due to the non-Abelian nature of Yang-Mills theories, it is highly non-trivial 
to verify diagrammatically the transversality of the gluon propagator, namely q v Yl llv {q) = 0. The 
main reason for this resides in the nature of the Slavnov-Taylor identities satisfied by the vertices 
appearing in the SDEs. These always involve auxiliary ghost Green's functions, and are either 
extremely complicated (four-gluon vertex) or cannot be cast in a useful form (ghost-gluon vertex). 
As a result, it is hard to determine which set of diagrams can be safely discarded in the gluon 
SDE (or in other words how to truncate the SDE) without violating the transversality property (i.e., 
the gauge invariance) of the answer. Keeping only diagrams (ai) and (02) is not correct even at 
one-loop, and adding (03) will not improve the situation beyond one-loop. 

To overcome this limitation, one needs to devise a truncation scheme for the non-Abelian 
SDEs that respects gauge invariance at every level of the dressed-loop expansion. This has been 
accomplished in [J3j] where the pinch technique has been used to carry out a systematic rearrange- 
ment of the entire Schwinger-Dyson series. For the case of the gluon self-energy, this procedure 
results in a new SDE describing a modified Green's function, namely the PT-BFM self-energy 
Yl^ v (q). As shown in Fig. ^ this new SDE still contains the conventional gluon propagator but is 
composed by new vertices (indicated with a tilde in what follows) which corresponds precisely to 
the one to be found when quantizing the theory within the BFM: the external gluons have been pro- 
jected dynamically into background gluons. This projection has far reaching consequences: when 
hit with the physical momentum q the modified (BFM) vertices will now satisfy Abelian Ward 
identities, implying in turn that one- and two-loop gluon and ghost contributions are individually 
transverse [Q, 

Notice that we do not have yet a dynamical equation for the conventional gluon propagator, 
since the quantity that appears at the lhs is not the one propagating in the diagrams of the rhs. 
This is solved by resorting to powerful identities known under the name of background quantum 
identities p2[], which relates rc-point Green's functions with different content of background and 
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Figure 4: The SDE for the gluon self-energy in the PT-BFM framework. 



quantum fields. In the case at hands one has 

A-V) = [1+G(<7 2 )] 2 A-V), (3-D 

where A(q 2 ) is the background gluon propagator and G(q 2 ) is the g^ v form factor appearing in the 
Lorentz decomposition of the auxiliary Green's function A^ v (q) (see Fig. |5|) 

A^v(^) = g 2 C A J D(k + q)A°(k)H av (-k,-q,k,q) 

= g ^ vG (q 2 ) + q -^L(q 2 ). (3.2) 
q 

The function H^ v (Fig. |5] again) is in fact a familiar object, since it appears in the all-order Slavnov- 
Taylor identity satisfied by the standard three-gluon vertex. It is also related to the full gluon-ghost 
vertex by the identity q v H^ v (—k — q,k,q) = — iT jU (— k — q,k,q); at tree-level, Hny = ig av and 
(~k,-q, k, q) = F fl (-k-q,k,q) = -q u . 

Considering only one-loop dressed diagrams (that is the only first two blocks in Fig. |J) one 
has then the following system of integral equations 4 

A-V) = * + ^ b f , 
iD-\p 2 ) = p 2 + iA Jv a {-k- p,p,k)A^ v (k)r v (-k 1 -p,k + p)D(p + k), 
J AP a {k)H ap (-k-q,k,q)D(k + q) 



G(q 2 ) ^ 



d-l 
1 

q 



+ i-j fq p A p(7 (k)r a (-k-q,k,q)D(k + q) 
q L Jk 



(3.3) 



One should be aware of the fact that there is no a-priori guarantee that the one-loop dressed gauge-invariant subset 
kept capture necessarily most of the dynamics, or, in other words, that they represent the numerically dominant contri- 
butions (however, for a variety of cases it seems to be true, see next section). But, the point is that one can systematically 
improve the picture by including more terms, without worrying that the initial approximation is plagued with artifacts, 
originating from the violation of gauge invariance. 
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Figure 5: Diagrammatic representation of the functions A and H appearing in the PT-BFM framework. 



where in the first term the transverse projectors have been traced out by virtue of the block-wise 



transversality of the PT-BFM equation, and the identity ( p.l\ ) have been used, to trade the BFM 
propagator for the conventional one. 

Before attempting to solve the system, there is one last step to perform, that is the projection to 
the Landau gauge (since we want to compare the SDE results to lattice simulations). This is a subtle 
exercise in the one-loop dressed gluon sector, for one cannot set directly = in the integrals due 
to the terms proportional to ^ l appearing in the BFM three-gluon vertex. Instead, one has to use 
the expressions for general E, , carry out explicitly the set of cancellations produced when the terms 
proportional to % generated by the identity k^A pv (k) = —i^k v /k 2 are used to cancel l/£ terms, 
and set <^ = only at the very end. It is relatively easy to establish that only the bare part of the full 
vertex contains terms that diverge as E, — > ||5||. One has 

= d^lUl r f (qA-k-q)A l ap (k)Al (J (k + q)L^( q ,k-k-q) 

+ I A PM [r^ ff (q, k, -k - q) + L»P° (q, k, -k - q)} - ( d - 1 f j\{k) 



(k + q) 7 
k»(k + q)n 



(3.4) 



i k k 2 {k+ q y y 

where L satisfies the Ward identity 

q^L^ pa (q, k, —k — q) = A _1 (& + q)P pa [k + q) - A" 1 (k)P pa (k) . (3.5) 
For the ghost diagrams one finds finally 

= ^ J^(~ k ~ 9' 9' k)D(k)D(k + qftni-k, -q, k + q)- 2idX jr>{k) . (3.6) 

4. Approximations, numerical results and comparison with lattice data 

Before attempting to solve the PT-BFM equations derived in the previous section, one needs 
to identify suitable approximations for the various vertices appearing in them. As far as the con- 
ventional gluon-ghost vertex entering in the ghost equation, lattice studies (in d = 4) indicate that 
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it deviates very mildly from its tree-level value which will be therefore kept throughout (the same 
tree-level approximation will be used for the auxiliary vertex function H av ). Instead, the strategy 
one should adopt for the background three-gluon and ghost-gluon vertices entering into the gluon 
propagator SDE, is the following. Given that the proposed truncation scheme relies crucially on 
the validity of the Ward identities satisfied by the background vertices, one should start out with 
an approximation that manifestly preserve them. The way to enforce this, familiar to the SDE 
practitioners already from the QED era, is to resort to the gauge-technique [[23|], namely solve the 
Ward identities. Specifically, one must express the vertices as a functional of the corresponding 
self-energies, in such a way that (by construction) the corresponding Ward identity is automati- 
cally satisfied 5 . One possibility, that contains the massless poles needed for the triggering of the 
Schwinger mechanism is [||] 

L p , pa (q,k,-k-q) = Y iipa (q,k,-k-q) +i% [U p(J (k + q) -U pa {k)] , 

(q, k, ~k-q)= f „ (q, k-k-q)-x% \t(k + q)- l(k)] . (4. 1) 

q 

Other forms are also allowed; in particular in a vertex that avoids the seagull divergences that 



the expressions ( |4.1| ) are bound to generate has been derived. 

The methodology described above constitutes, in fact, the standard procedure even in the con- 
text of QED, where the structure of the SDE is much simpler (the SDE for the photon contains one 
single graph that involves the photon-electron vertex which satisfies automatically a Ward identity). 
Thus, while the PT-BFM approach described here replicates QED-like properties at the level of the 
SDEs of QCD, which is a striking fact in itself, does not make QCD easier to solve than QED. One 
should appreciate an additional point though: any attempt to apply the approach described above 
in the context of the conventional SDEs is bound to lead to the violation of the transversality of 
the gluon self-energy, because (i) the vertices satisfy non-linear Slavnov-Taylor identities (a fact 
that makes the application of the gauge-technique impractical), and (ii) even if one had managed 
to come up with good Ansatze for all vertices, one should still keep all self-energy diagrams in 
Fig. H to guarantee that q v Yl^ v (q) = 0. From this point of view, the improvement of the PT-BFM 
approach over the standard formulation becomes evident. 

4.1 Gluon and ghost Green's functions 



The solutions of the PT-BFM equations (3.3) with the vertices (4.1) in the SU(3) d = 4 



case [24], and in the one-loop approximation for the SU(2) d = 3 case [25], are shown in Fig. 
As it can be seen in the plots, the PT-BFM equations capture qualitatively the correct d = 4 physics 
even within the rough approximations used, while in the d = 3 case the agreement becomes quan- 
titative. In this latter case one gets the value m/2g 2 = 0.153 in agreement with many independent 
lattice studies. 

The agreement found between the SDEs and the lattice results allows one to study other quan- 
tities of interest by using the lattice directly as an input into the various SDE. The general strategy 



5 Notice that this procedure leaves the transverse (i.e., automatically conserved) part of the vertex undetermined. 
This is where the SDE for the vertex enters: It is used precisely to determine these transverse parts. Since however the 
transversality of the self-energy depends on the longitudinal part of the vertex only, one can apply all sort of approxima- 
tions for the determination of the transverse form factors. 
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Figure 6: Gluon propagator and ghost dressing function for the SU(3) d = 4 and SU(2) d = 3 cases, and 
comparison with the corresponding lattice data. 



adopted in this case is the following. One takes the lattice gluon propagator as an input for the 
ghost SDE; then solves for the ghost dressing function, tuning the coupling constant g in such a 
way that the solution gives the best possible approximation to the lattice result. Obviously one must 
check that the coupling so obtained (at the renormalization scale used for the computation) is fully 
consistent with known perturbative results (obtained in the MOM scheme, which is the scheme 
used in our computations). Then the system is "tuned", and one can construct and analyze other 
quantities built out of A, F and g, as is done in the next two sections. 

4.2 The G and L auxiliary functions 

The first quantities that can be studied by means of the procedure just described, are the aux- 
iliary functions G(q 2 ) and L(q r), which, within our approximations, are given by 



G(q 2 ) 
L(q 2 ) 



g 2 c A 

d-ljk 
8 2 Ca f 
d-lJk 



(d-2) + 



(k-qf 



lc 2 q 2 



A(k)D(k + q), 



l-d 



(k-q) 7 
k 2 q 2 



A(k)D(k + q). 



(4.2) 



Before doing that, let us notice that in the Landau gauge one can prove that the dressing function 
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Figure 7: The form factors —G(q 2 ) and L(q 2 ) determined from Eq. (4.2) at different normalization points 
fl, using the procedure specified in the text. In the bottom panels we show the direct calculation on lattice of 
the Kugo-Ojima function u{q 2 ) at fi = 4 GeV and a comparison with our result. 



F{q ) and the form factor G(q ) and L(q ) are related through the BRST identity [26] 

F- l (q 2 ) = \+G(q 2 )+L(q 2 ). (4.3) 

Since, under very general conditions on the gluon and ghost propagators, L(q 2 ) — > when q 2 — > 
one has the IR relation F (0) = 1 +G(0). Thus, we see that a divergent (or enhanced) dressing 
function requires the condition G(0) = —1. The latter looks suspiciously similar to the Kugo- 
Ojima confinement criterion demanding that a certain function u{q 2 ) (the Kugo-Ojima function) 
acquires the IR value w(0) = — 1. Indeed it is possible to show that G is nothing but the Kugo-Ojima 



function [26] 



u{q l )=G{q l ), (4.4) 

and therefore G{q 2 ) encodes practically all relevant information of the IR dynamics of the ghost 
sector. 

In Fig. ([7j) we show the auxiliary functions (4.2) calculated in d = 4 at different renormalization 



points. One can see that indeed L(0) = and that in general L(q 2 ) is suppressed with respect to 
G{q 2 ) (also, power counting shows that L is UV finite, while G is not). In addition one finds that 
the function G(q 2 ) saturates at an IR value bigger than —1 (around —2/3 for the renormalization 
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points chosen) once again excluding IR enhancement of the ghost dressing function. These results 
can be compared with the direct lattice calculations of the Kugo-Ojima function in terms of Monte 
Carlo averages [|27|], though the comparison is at most suggestive because of systematic errors due 
to the different renormalization procedures employed. 

4.3 The effective charge 

Another important information that can be extracted from the PT-BFM equations is the run- 
ning of the renormalization group (RG) invariant QCD effective charge, and, in particular, its be- 
havior and value in the deep IR region. This quantity lies at the interface between perturbative 
and non-perturbative effects in QCD, providing a continuous interpolation between two physically 
distinct regimes: the deep UV, where perturbation theory is reliable, and the deep IR, where non- 
perturbative techniques must be employed. 

There are two possible RG invariant products on which a definition of the effective charge 
can be based: r~(q 2 ) which exploits the non-renormalization property of the ghost vertex in the 
Landau gauge, and d(q 2 ) which exploits the fact that BFM quantities satisfy Ward (as opposed to 
Slavnov-Taylor) identities. One has 

?{q 2 ) = *V)A(«rW) d(q 2 ) = g 2 (^L 2 )A(q 2 ). (4.5) 

These two dimensionful quantities, that have a mass dimension of —2, share an important common 
ingredient, namely the scalar cofactor of the gluon propagator A(q 2 ) which actually sets the scale. 
The next step is then to extract a dimensionless quantity that would correspond to the nonpertur- 
bative effective charge. Perturbatively, i.e., for asymptotically large momenta, it is clear that the 
mass scale is saturated simply by q 2 , the bare gluon propagator, and the effective charge is defined 
by pulling a q~ 2 out of the corresponding RG-invariant quantity. Of course, in the IR the gluon 
propagator becomes effectively massive; therefore, particular care is needed in deciding exactly 
what combination of mass scales ought to be pulled out. The correct procedure in such a case [||] 
is to pull out a massive propagator of the form (in Euclidean space) [q 2 + m 2 (q 2 )]~ l , with m 2 {q 2 ) 



the dynamical gluon mass 6 . One then has, making use of the background quantum identity (3.1) 

a gh (q 2 ) = a(LL 2 )[q 2 +m 2 (q 2 )]A(q 2 )F 2 (q 2 ); a(q 2 ) = a^ 2 )[q 2 +m 2 (q 2 )] ff^ 2)]2 - (4- 



6) 



In addition, due to the identity (|4.3|), the two effective charges are related through [ |13| ] 

L(q 



a(q l ) = a gh (q z 



^ n2 



1 + 



l+G(q 2 ) 



(4.7) 



Since L(0) = 0, we therefore see that not only the two effective charges coincide in the UV region 
where they should reproduce the perturbative result, but also in the deep IR where one has a(0) = 
a g h(0). In addition, due to the relative suppression of the L{q 2 ) form factor as compared to G(q 2 ) 
even in the region of intermediate momenta (Fig. [7]), where the difference reaches its maximum, 

6 Given that the gluon propagator is finite in the IR, if one insists on factoring out a simple q~ 2 term, one would 
get a completely unphysical coupling, namely, one that vanishes in the deep IR, where QCD is expected to be (and is) 
strongly coupled. 
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Figure 8: Left panel: Comparison between the two RG-invariant products d(q 2 ) (solid line) and ?(g 2 ) 
(dashed line); notice that there are two overlapping curves at different fi for each product. Right panel: 
Comparison between the QCD effective charge extracted from lattice data: cc(q 2 ) (red line with circles) and 
0!gh (black line with squares) for two different masses: mo — 500 MeV (dashed) and mo — 600 MeV (solid). 



the relative difference between the two charges is small (at most 10%), making them practically 
indistinguishable. 

In Fig. H we show both a check of the RG-invariance of the combinations (03) as well as 



a comparison between the effective charges (i.t) extracted from the lattice data for two different 
values of the running gluon mass [13]). 



5. Conclusions and outlook 



In this talk we have reviewed the PT-BFM framework for formulating (and solving) the SDEs 
of Yang-Mills theories. In particular, we have sketched how a dynamically generated gluon mass 
gives rise self-consistently to condensation of (thick) vortices leading to confinement, and thus that, 
contrary to naive expectations, a gluon mass does not conflict with confinement, but enables it. In 
addition, we have derived within the PT-BFM framework a new set of SDEs that (i) accommodate 
the dynamical gluon mass generation through the Schwinger mechanism, and (ii) have much better 
truncation properties than the conventional equations. The solutions of this new set of SDEs for 
the cases of the gluon propagator and the ghost dressing function shows IR saturation in agreement 
with all lattice simulations up to date. The analysis of the effective charge has also revealed that the 
gluon mass keeps the theory strongly coupled and allows for the presence of a conformal window. 

It should be emphasized that within this framework, and contrary to other approaches, the 
BRST symmetry of the theory remains intact throughout. Indeed, at no point have we tampered 
with the original Yang-Mills Lagrangian (and therefore we have no tree-level BRST violation), 
while the special transversality properties of the new SDEs ensure the preservation of the STIs (a 
direct consequence of the BRST symmetry), and most importantly of the transversality condition 
q v Yl llv = 0, at every step of the truncation. 

There are two complementary aspects that need to be further investigated. On the one hand, 
one would like to improve the agreement between the PT-BFM and lattice results in d = 4. This is in 
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our opinion linked to devising better vertex Ansatze for implementing the Schwinger mechanism. 
On the other hand, the lattice should address the issue of how the results changes when calculations 
are performed in gauges other than the Landau (recently an algorithm implementing R^-type of 
gauges on the lattice has been developed, see [^]). In particular, a lattice implementation (along 



the lines suggested in [29] of the BFM in the Feynman gauge, where all results are known to be 



free from gauge artifacts [§, 16], 17], would be more than welcome. 

Finally, it is now timely to focus also on phenomenological issues in order to complement the 
lattice findings with experimental evidence. For example in [ ]30[ ] the influence of an IR dynamical 
gluon mass scale in the calculation of pp and pp forward scattering quantities through a QCD- 
inspired eikonal model was addressed, finding good agreement with data for a typical gluon mass 
of the order of 500 MeV. Also for hybrids with an extra gluon it was shown in [ pT| ] that the best fit 
to lattice data on three heavy quarks plus a gluon definitely favor a gluon mass of the same order. 

It is evidently this interplay between theory, lattice simulations and experiments that in the 
long run will prove to be most fruitful in unraveling the full IR dynamics of Yang-Mills theories. 

Acknowledgments: We would like to thank the organizers of LC2010 for providing such a stimu- 
lating environment and enjoyable workshop. 
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